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Abstract 


2-  /?  /  ^ b,  r 

We-consider  the  random-accessing  problem  of  a  single. 


collis ion- type ,  slotted,  packet-switched  communication  channel 


by  a  large  number  of  independent,  data  transmitting  bursty 
users.  |i«r propose  and  analyze  an  easy-to-implement  algorithm 


under  the  realistic  assumption  that  each  user  inspects  the 
channel  outcome  feedback  only  whenever  he  is  blocked.  Wa-  aooumo  • 


binary  feedback  which  informs  the  users  only  about  whether  or  not 

jLy  \  Ao  rs* 

there  was  a  collision  in  the  previous  slot.  Me  a  hew--  that  the 


algorithm  results  in  finite  average  delays  for  transmission  at 
rates  less  than  0.36  packets  per  channel  slot,  and^we  give  an 
exact  upper  bound  for  the  average  delay/  ;s  a/vir. 


1.  Introduction 


The  multiple-access  problem  in  communications  is  the  problem 
of  organization,  or  coordination  of  a  population  of  users  for 
the  efficient  sharing  of  the  resources  of  a  single  channel  used 
by  the  users  for  information  transmission.  This  situation 
arises  in  a  number  of  applications :  Computer-communication 
networks,  packet-radio  networks,  satellite  communication  networks, 
local  area  networks . 

Random  multiple-access  schemes  are  an  important  class  of 
techniques  that  employ  distributed  control  algorithms  to  cope 
with  the  multiple  access  problem.  These  schemes  are  especially 
useful  in  the  presence  of  an  asymptotically  large  number  of 
ill-specified, independent  bursty  users.  The  users  gain  access 
into  the  channel  on  a  contention  basis.  The  accessed  channel  is 
a  collision- type,  packet  switched,  time-slotted  transmission 
channel.  Some  form  of  feedback  information  associated  with 
the  message  transmissions  is  always  assumed  to  be  available  to 
the  contenting  users.  For  this  general  model,  a  variety  of 
access  algorithms  (protocols)  has  been  proposed  and  analyzed 
by  several  authors.  The  properties  of  these  algorithms  vary 
considerably  with  the  level  of  the  feedback  information  assumed 
available  to  the  users. 

In  [4-8]  it  is  assumed  that  immediately  after  each  channel 
slot  a  ternary  feedback  is  broadcasted  to  the  users.  This  is 


known  as  0,l,e  feedback  and  informs  the  users  whether  the  previous 
slot  was  empty  (0)  ,  or  contained  one  packet  (1) ,  or  contained  a 
collision  (e) .  A  collision  occurs  wnenever  more  than  one  users 
attempt  transmission  within  the  same  slot.  All  information 
contained  in  the  collided  packets  is  assumed  lost,  and  these 
packets  must  be  retransmitted  at  later  times. 

The  algorithms  developed  in  [9]  use  binary  feedback.  Binary 
feedback  is  less  informative  compared  to  ternary  feedback  and 
may  be  available  in  three  different  forms:  "Collision/No  Collision" 
feedback,  "Something/Nothing"  feedback  and  "Success/Failure" 
feedback  (notation  suggested  by  Mehravari  and  Berger  [9]). 

Several  recent  efforts  of  developing  more  efficient 
realizable  algorithms  have  used  more  informative  types  of  feedback 
than  ternary  feedback.  In  {10]  it  is  assumed  that  after  each 
collision  the  number  (up  to  an  upper  maximum  limit)  of  the 
packets  involved  is  revealed  to  all  users,  through  a  bank  of 
energy  detectors.  Also,  in  111]  it  is  assumed  that  additional 
information  (four-valued,  or  five-valued  feedback)  is  available 
to  the  users  through  the  use  of  control  mini-slots.  All  the 
algorithms  in  the  papers  mentioned  so  far  require  that  each  user 
inspects  the  feedback  broadcasting  for  every  channel  slot  over 
the  entire  operation  of  the  random-access  system.  In  the  slotted 
Aloha  algorithm  [1,  2]  (which  is  unstable  for  an  asymptotically 
large  number  of  users)  each  user  inspects  only  the  slots  that 
correspond  to  his  own  attempts. 

Tsybakov  and  Vvedenskaya  [12]  proposed  and  analyzed  a 
"limited-sensing"  algorithm,  called  "Stack"  algorithm,  where 
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each  user  inspects  the  feedback  broadcasting  only  whenever  there 
is  an  unsuccessfully  transmitted  packet  in  his  buffer.  The 
"Stack"  algorithm  uses  ternary  feedback  and  achieves  a  maximum 
stable  throughput  of  at  least  0.384  packets  per  slot. 

Using  the  same  feedback  level  as  in  [12]  Papantoni-Kazakos 
and  Marcus  [13]  developed  a  limited  channel  sensing  algorithm 
for  a  limited  number  of  data  users. 

In  the  present  paper,  we  propose  and  analyze  a  Limited 
Sensing  random  access  Algorithm  with  Binary  Feedback  (LSBFA) . 

The  user  and  channel  models  assumed  are  described  in  section  2. 

In  the  LSBFA,  users  with  new  packets  transmit  their  packets  in 
the  first  slot  following  their  arrival,  and  then  they  resolve 
any  collisions  using  the  "Capetanakis-Tsybakov-Mikhailov- 
Collision-Resolution-Algorithm"  (CTMCRA)  [14].  In  contrast  to 
the  "Tree- type"  random  access  algorithms  [4-11] ,  where  there  is 
an  explicit  and  separate  collision  resolution  period,  the 
LSBFA,  like  the  "Stack"  algorithm,  allows  new  packets  to  continuously 
enter  into  the  system  independently  of  the  collision  resolution 
process  already  in  progress.  From  a  practical  point  of  view,  the 
"continuous-entry"  feature  is  very  significant,  since  the  users 
monitor  the  feedback  channel  only  whenever  they  have  a  packet  to 
send  (limited-sensing).  Thus,  the  undesirable  in  several 
applications  necessity  of  all  users  monitoring  the  feedback 
channel  constantly — even  if  they  have  no  packet  to  send — is 
eliminated. 

The  organization  of  the  paper  is  as  follows: 


In  section  2,  we  present  the  user  and  channel  models,  and 


the  LSBFA  statement  and  general  operation.  j 

j 

In  section  3,  we  analyze  the  algorithm  by  expressing  and 
studying  a  system  of  recursive  equations  for  the  expected 
length  of  collision  resolution  sessions;  we  evaluate  the 
stability  region  of  the  algorithm  and  the  expected  length  of 
a  session. 

In  section  4,  we  use  the  results  of  section  3  to  give 
an  exact  upper  bound  for  the  average  packet  delay. 

In  section  5 ,  we  compare  the  LSBFA  to  other  random  access 
algorithms. 


2.  The  Model  and  the  Statement  of  the  Algorithm 

We  assume  that  an  infinite  population  of  spatially 
isolated,  independent,  bursty  users  (transmitters)  share  a 
single  channel  to  communicate  with  a  central  facility  (common 
receiver) .  The  users  transmit  data  packets  of  fixed  duration 
taken  to  be  the  unit  of  time.  The  channel  time  is  divided 
into  unit-time-segments  called  slots.  The  unit  interval 
(t,  t+1)  is  called  slot  t  (t  =  0 ,  1 ,  2 ,  .  .  . ) .  All  users 
are  synchronized  to  the  starting  points  of  the  channel  slots, 
and  they  attempt  transmission  of  some  packet  only  at  the 
beginning  of  some  channel  slot. 

The  users  do  not  communicate  with  each  other  directly; 
therefore,  the  probability  of  more  than  one  users  attempting 
transmission  of  a  packet  within  the  same  slot  is  nonzero.  A 
channel  slot  is  a  collision  slot  if  more  than  one  packets 
attempted  transmission  within  it.  All  information  in  the 
packets  involved  in  a  collision  is  assumed  lost,  and  these 
packets  have  to  be  retransmitted.  A  channel  slot  is  empty  if 
no  packet  attempted  transmission  within  it,  and  it  is  successful 
if  exactly  one  packet  was  transmitted  within  it.  In  the  later 
case  it  is  assumed  that  the  transmitted  packet  reaches  its 
destination  error  free. 

A  feedback  channel  from  the  common  receiver  informs  the 
transmitters  at  the  end  of  each  slot  whether  or  not  there  was 


a  collision  in  that  slot.  We  assume  that  the  feedback  channel 


is  a  noiseless  broadcast  channel,  and  that  propagation  delays 

are  negligible.  Therefore,  immediately  at  the  end  of  slot  t 

a  user  who  is  interested  in  the  outcome  of  that  slot  can  learn 

it  be  monitoring  the  feedback  channel.  Let  the  random 

variable  Z.  denote  the  outcome  of  slot  t.  We  have 
t 


INC  if  slot  t  contained  £  1  packets  (no  collision) 
C  if  slot  t  contained  >_  2  packets  (collision) 


The  Collision/No  Collision  (CNC)  binary  feedback  assumed 
here  uses  less  information  compared  to  the  0,1, e  ternary 
feedback,  since  the  former  does  not  distinguish  an  empty  slot 
from  a  successful  slot.  Its  implementation  can  be  based  on 
a  simple  binary  acknowledgment  scheme  from  the  central  facility 
to  the  users  ("NC"  or  "C") . 

Let  the  random  variable  N(t)  denote  the  number  of  new 
packets  appearing  in  the  system  for  transmission  from  all  users 
combined  during  slot  t.  It  is  assumed  that  {N(t)}(t  =  0,  1,  2,  .  .  . 
is  a  sequence  of  independent  and  identically  distributed 
random  variables.  Let  pn  =  Pr{N(t)  =  n}  be  the  probability 
mass  function  (p.m.f.)  of  N(t)(t  =  0,  1,  2,  .  .  .),  and  let 
X  =  E (N ( t) }  be  its  expectation,  which  is  assumed  being  finite. 

Thus,  X  is  the  intensity  of  the  cumulative  input  traffic 
measured  in  number  of  packets  per  slot.  For  the  infinite-population 
model,  only  one  packet  requiring  transmission  can  be  present 
at  a  station  at  any  given  point  of  time. 

The  following  definition  will  be  used  in  the  statement  of 


the  algorithm. 
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Definition  1  At  any  time  t  a  user  may  be  either  active  or  inactive 
At  any  time  t  an  active  user  may  be  either  new 
or  blocked 

A  new  user  at  time  t  is  a  user  with  a  packet 
generated  during  slot  t-1 

A  blocked  user  at  time  t  is  a  user  with  a  packet 
that  has  attempted  transmission  and  experienced 
collision  at  some  slot  prior  to  slot  t 
At  time  t  a  packet  is  new  or  blocked  if  it  belongs 
to  a  new  or  blocked  user  respectively 

The  statement  of  the  algorithm  is  contained  in  two  simple 
rules.  For  the  implementation  of  the  algorithm  in  a  distributed 
fashion,  it  suffices  for  each  user  to  have  a  counter  and  a 
binary  fair  coin.  The  rules  of  the  algorithm  are  followed  by 
the  active  users  only,  and  are  as  follows : 

Let  the  random-access  system  start  at  t=0  with  all  counters 
set  at  0. 

Rule  1  At  time  t  (t=0,  1,  2  .  .  .)  blocked  users  with 
counter  at  0  and  all  new  users  transmit  their 
packets  in  slot  t 

Rule  2  At  the  end  of  slot  t  (just  prior  to  time  t+1)  all 
active  users  inspect  the  feedback  channel.  If 
slot  t  were  collison  free  (Zfc  =  NC) ,  the  user  who 
transmitted  his  packet  (if  any)  leaves  the  system 
(becomes  inactive) ,  and  all  blocked  users  decrement 
their  counters  by  one.  If  slot  t  were  a  collision 
slot  (Zfc  =  C) ,  each  collided  user  tosses  a  binary 


1 


fair  coin  and  sets  his  counter  to  0  or  to  1 
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according  to  the  outcome  of  the  coin  tossing.  All 
other  blocked  users  increment  their  counters  by 
one . 


i 


] 


3 .  Algorithm  Analysis 


The  two  most  important  performance  measures  of  a  random- 
access  algorithm  are  the  average  delay  of  a  packet  and  the 
maximum  stable  throughput  of  the  system.  The  delay  of  a  packet 
is  the  time  between  the  instant  the  packet  originates  as  a 
new  packet  until  the  instant  it  is  successfully  transmitted. 

til 

Let  6n  be  the  random  variable  that  denotes  the  delay  of  the  n*" 
packet. 

A  random-access  algorithm  or  system  is  called  stable  if 
the  lim  supE{6n)  is  finite,  assuming  that  the  limit  exists. 

This  means  that  for  a  stable  algorithm  the  delay  of  a  packet 
will  remain  finite  with  probability  one.  The  throughput  or 
output  rate  of  a  random-access  system  is  the  long-run  average 
number  of  successfully  transmitted  packets  per  unit  time;  it 
is  denoted  by  n-  Given  an  algorithm,  let  n*  be  the  supremum 
of  n-  If  the  input  rate  X  is  less  than  n*,  then  the  throughput 
is  x  and  the  system  is  stable.  If  the  input  rate  x  exceeds  the 
maximum  output  rate  n*,  then  the  average  packet  delay  becomes 
unbounded  (Little's  result)  and  the  system  is  unstable.  Thus, 

n*  =  sup{X:  the  system  is  stable) 

We  call  n*  the  efficiency  of  the  given  algorithm.  The 
interval  (0,  n*)  is  its  stability  region. 
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In  this  section  we  study  the  stability  region  and  the 
average  packet  delay  for  the  algorithm  described  in  section  2. 

We  proceed  with  the  following  definitions: 

Definition  2  A  clock  instant  tp  is  a  renewal  instant 

if  all  active  users  at  tR  (if  any)  are  new 
users 

Definition  3  A  session  is  a  sequence  of  consecutive 

channel  slots,  that  begins  and  ends  at  two 
consecutive  renewal  instants. 

Definition  4  A  session  starting  at  some  renewal  instant 

t„  is  called  a  session  of  multiplicity  k 
K  " 

if  the  number  of  active  (new)  users  at  tp 
is  k;  k  =  0,  1,  2,  .  .  . 

If  t_  is  a  renewal  instant,  then  at  tD,  by  definition  2, 
there  are  no  blocked  users  in  the  system.  Thus,  any  previously 
blocked  packet  has  been  successfully  transmitted  prior  to  slot 
tR.  From  definitions  1  and  3,  and  the  independence  of  the 
incoming  traffic  from  a  particular  session,  it  is  clear  that  the 
lengths  of  consecutive  sessions  are  independent  and  identically 
distributed  random  variables.  Let  the  random  variable  t  denote 
the  length  of  an  arbitrary  session,  and  let  denote  the  length 
of  a  session  of  multiplicity  k.  The  distribution  of  both  random 
variables  t  and  depends  only  on  the  probability  mass  function, 
that  models  the  incoming  traffic  (input) ,  and  on  the  rules  of  the 
algorithm,  but  not  on  the  particular  session. 

The  sessions  with  multiplicity  0  or  1  are  trivial  and  both 
have  length  equal  to  one  slot.  The  session  of  multiplicity  0  is 
called  the  empty  session. 
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Let  tg  and  t  be  the  two  random  consecutive  renewal  instants 
that  denote  the  starting  and  ending  instants  of  the  nfc^  nonempty 
session  respectively.  Since  at  both  t  and  t  there  are  no 
blocked  users  present  at  the  system,  the  number  of  successfully 
transmitted  packets  during  the  course  of  the  session  is  simply 
the  random  number  of  packets  that  appeared  to  the  system 
requiring  transmission  during  the  time  interval  (t  -1,  t  -1) . 

Let  the  random  variable  M  denote  the  random  number  of  success¬ 
fully  transmitted  packets  during  the  nth  nonempty  session.  The 
average  packet  delay  of  a  nonempty  session  is  defined  as  follows: 

1  M 

D  =  E{M  l  6)  (1) 

-ill  m 

m=l 

where  6^  is  the  delay  experienced  by  the  mth  successfully 
transmitted  packet  during  the  session. 

It  is  clear  that  D  is  independent  of  the  particular  session, 
because  sessions  are  independent  of  each  other,  and  because  the 
input  traffic  is  a  process  with  independent  and  identically 
distributed  increments  (N(t)5.  Thus,  the  average  packet  delay 

4* 

of  the  n  session,  given  by  (1) ,  is  the  average  packet  delay 
of  the  system,  that  is  over  all  sessions. 

Since  all  the  packets  of  the  session  requested  transmission 
not  earlier  than  t  ,  and  were  successfully  transmitted  not  later 
than  t  -1,  we  have 

<  t-1;  m  =  1,  2,  .  .  .  (2) 

m  — 

where  t  is  the  length  of  the  session.  From  (1)  and  (2)  we  have 


the  following  upper  bound  for  D: 
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The  random  variable  I  is  binomially  distributed: 

Pr{ I=i }  =  bk(l)  =  (k)2_k 

The  random  variables  M  and  N  are  identically  distributed  with 
Pr{M=i}  =  Pr{N=i}  =  p(i) 

where  p(*)  is  the  p.m.f.  of  the  input  increment. 

Proof 

For  the  trivial  sessions  of  multiplicity  0  or  1,  by  definition, 

Tg  =  t.  =  1.  For  k  _>  2  let  the  session  start  at  the  renewal 

instant  ts  with  k  new  users.  According  to  the  first  rule  of  the 

algorithm  all  the  k  users  transmit  their  packets  in  slot  tg. 

Hence,  slot  tg  is  always  a  collision  slot.  At  instant  tg+l  there 

are  no  other  blocked  users  except  for  the  k  collided  users,  who, 

according  to  the  second  rule,  toss  a  fair  coin  and  set  their 

counters  either  to  0  or  to  1.  Let  the  random  variable  I  denote 

the  number  of  those  from  the  k  users,  who  set  their  counters  to 

0.  Then,  k-I  is  the  random  number  of  blocked  users  with  their 

counters  set  at  1.  Clearly,  Pr{I=i}  =  (^)2  ;  i  =  0,  1,  .  .  ,,k. 

Let  the  random  variable  M  denote  the  number  of  new  users 

at  instant  t  +1.  These  are  the  users  with  a  packet  originated 

during  slot  t  .  Clearly,  Pr{M=m}  =  Pr{N(t  )  =  m}  =  p  ,  and  M 
s  s  m 

is  independent  of  1.  According  to  the  algorithm  a  total  number  of 
I+M  users  transmit  their  packets  in  slot  tg+l,  while  the  k-I 
blocked  users  with  their  counters  previously  set  to  1  inspect  the 
feedback  channel  and  increment  (decrement)  their  counters  by  one 
for  each  subsequent  collision  (collision-free)  slot  respectively. 
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The  crucial  observation  here  is  that  all  the  k-I  blocked 

users  have  identical  counter  indication  until  the  random  instant 

t  at  which  this  indication  becomes  0  for  the  first  time, 
o 

Furthermore,  it  is  not  difficult  to  see  from  the  rules  of  the 
algorithm,  that  the  identical  counter  indication  of  the  k-I 
blocked  users  is  always  greater  than  the  counter  indication  of 
any  other  blocked  user  in  the  system.  Thus,  at  instant  t  there 
are  no  other  blocked  users  in  the  system  except  of  the  k-I 
users  who  all  have  their  counters  set  to  0.  This  means  that  the 
I+M  packets  that  started  accessing  the  channel  in  slot  tg+^ 
and  all  the  packets  that  appeared  to  the  system  during  (t  -1, 
t  -1)  have  been  successfully  transmitted  by  the  random  instant 

V 

Let  the  random  variable  N  denote  the  number  of  new  users 

at  t  .  Clearly,  Pr{N=n}  =  Pr{N(t  -1)  =  n}  =  p  and  N  is 
o  u  o  n 

independent  of  I  and  M.  In  slot  tQ  a  total  of  k-I+N  users 

transmit  their  packets,  since  there  are  k-I  blocked  users  with 

counter  indication  0  and  N  new  users.  Consequently,  the  I+M 

users  may  be  thought  of  as  starting  an  independent  session  of 

random  multiplicity  I+M  ,  that  begins  at  tg+l  and  ends  at  tQ. 

This  session  is  immediately  followed  by  a  session  of  multiplicity 

k-I+N.  This  later  session  starts  at  t  and  ends  at  t0,  which 

is  the  first  renewal  instant  after  t  .  Thus, 

s 


1  +  [to 


<v1)] 


(te  -  V 


=  1  +  T  +  T 

I+M  k-I+N 


Q.E.D. 
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Let  G{k,  z)  ,  0  <_  z  <_  1  denote  the  moment  generating  function 

of  the  random  variable  t. -Is 

k 

Tk-1 

G(k,  z)  =  E  (z  ) 

In  view  of  Theorem  1,  the  following  theorem  is  straightforward: 
Theorem  2 


G(0,  z) 
G (k ,  z) 


G  (1 ,  z)  *  1 

k  00  «  2 

III  b.  (i)  p  (m)  p  (n)  z  G(i+m,  z)G(k-i+n,  z) 
i=0  m=0  n=0 

(6) 


where  bk(i) 


(^)2  k  and  p( •) 


is  the  p.m.f.  of  the  input  increment. 


For  the  mean  length  of  a  session  of  multiplicity  k  we  have: 


Lk  '  E<lk>  -  Jz  G<k-  *> 


+1;  k  =  0 ,  1 ,  2 ,  .  .  . 


z=l 


By  differentiating  (6)  we  respect  to  z,  or  by  directly 
taking  expectations  in  (5)  we  have: 


Corollary  1 


k  k  °° 

=  1  +  2  I  I  b.(i)p(j)L.  .  =  1  +  2  l  q  (m)L  ;  k  >  2 

K  i=0  j=0  K  1  +  3  m=0  k  m 

(7) 

where 

v 

q^tm)  =  bk(m)*p(m)  =  l  bk(i)p(m-i)  (8) 

v  =  min(k,  m) 


♦denotes  convolution 
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It  is  noteworthy  that  the  coefficients  q^f*),  given  by  (8), 
are  the  values  of  the  convolution  of  the  binomial  p.m.f.  b^f*) 
with  the  p.m.f.  of  the  input  increments  p(*).  Thus,  q^f*)  is  the 
p.m.f.  of  the  sum  I+J  of  two  independent  random  variables:  The 
random  variable  I,  which  is  binomially  distributed,  and  the 
random  variable  J,  which  is  distributed  according  to  p(*). 

But  this  is  exactly  what  the  continuous  entry  algorithm  does 
to  resolve  collisions.  After  a  collision  the  collided  users  (k) 
toss  a  fair  coin, and  those  who  tossed  0  (I)  transmit  their 
packets  in  the  next  slot  along  with  the  newcomers  (J) .  Note 
also  that  for  p(0)  =  1,  p(i)  =  0,  i>0  equation  (7)  becomes 
equation  (3.12)  of  [14],  that  gives  the  system  of  equations 
satisfied  by  the  conditional  mean  length  of  a  collision  resolution 
interval  for  the  CTMCRA. 

The  system  of  linear  equations  for  L^,  given  by  (7) , 
will  be  of  central  interest  in  this  paper.  In  what  follows,  we 
investigate  the  conditions,  under  which  system  (7)  has  a  unique 
nonnegative  solution,  such  that  0  <  1^  <  "  for  0  _<  k  <  ®. 
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3. 2  Stability 

We  consider  a  general  system  in  k  that  corresponds  to 
the  system  given  by  (7) : 


1 


xk 


1  +  2  T  q,  (m)  x  ;  k  >  2 
m=0  *  m 


(9) 


where  (m)  are  as  given  in  (8). 

We  are  interested  in  investigating  the  conditions,  under 
which  system  (9)  has  a  nonnegative  solution,  that  is  bounded 
for  finite  k. 

Given  a  sequence  Y  =  {y^}  with  y^eR  (k  =  0,  1,  2,  .  .  .)» 
we  define  the  operators  A  and  B  (assuming  that  they  exist)  as 
follows : 


A [Y ]  =  { A^ [Y] } ;  k  *  0,  1,  2,  .  .  . 
such  as:  AQ [ Y I  =  A^fY]  =  1 

CO 

Akm  =1  +  2  l  <3k(ro>ym'*  k  1  2 
m=0 


and 


B [ Y ]  =  { Bk [Y] };  k  =  0,  1,  2,  .  .  . 
such  as:  Bg[Y]  =  B^[Y]  =  0 

Bk[Y}  =  2  l.  ^k(m)Ym;  k  -  2 

m=2 


(10) 


(ID 


Let  also  An [Y ] (Bn [Y ] )  denote  the  sequence  resulting  from 
the  n  times  repeated  application  of  operator  A(B)  respectively 
on  the  initial  sequence  Y  (A°(Y]  =  B°[Y]  =  Y) . 
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Using  the  operator  A  system  (9)  becomes: 


\  =  AkfX];  k  >  2 

The  following  theorem  gives  a  general  sufficient  condition, 
under  which  system  (9)  has  a  unique,  nonnegative,  bounded  (for 
finite  k)  solution. 

Theorem  3 

In  the  class  of  sequences  X  such  that 

oo 

lim  max  £  q,  (m)  |x  |  =  0  (12) 

i-*~  k<i  m=i  K  m 

system  (9)  has  a  unique  solution  X  =  (xk>  if  the  p.m.f.  of  the 
input  increment  p(-)  is  such  that,  there  exists  some  nQ  <  °° 

(n^  =  1,  2,  .  .  .),  such  that 

n0 

(B  )k[F(A,  r)  ]  >  0  for  every  k  2  (13) 

where 

F  ( A  ,  r)  =  ( 2_k  f  k  ( A  ,  r)  };  k  >  2 

f  k  ( A  ,  r)  =  (1  -  2  A  )  k  +  r  (1  -  2A)  -  2A  (14) 

OO 

A  =  £  ip(i),  r  =  p(l)/p(0),  and  B  is  the 

i=0 

operator  defined  in  (11) 

If  (13)  is  true,  then  for  every  k,  we  have: 

0  (An)k[*°^X]  xk  for  every  n  >_  0; 
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<_  (An)kfx^]  for  every  n  _>  nQ,  and 
\  =  lim{An)kI (0)X]  =  lim(An)k[X(0) ] 


where 


(0>X  =  {(0)  k}  with  <0)x0  =  (0) 


x^  =  1  and 


(0) 


xk  =  b 'k-c ' ; 


k  >  2 

b'  =  2/(l-2A)  ,  c'  =  2b'  A+l 

X(0)  =  {xk(0>>  with  Xq(0)  =  xj0)  =  1  and  x^0*  =  bk-c;  k  >  2 

b  =  max ( ( B  °)  [G(r)]/(B  °).[F(A,  r)  ] ) 

k>2  k 

G  (r)  =  (2~kgk(r)},  gk(r)  =  2(k+r+l) 

F (x ,  r)  is  as  defined  in  (14),  A  is  as  defined  in  (10) 
and  B  i«5  as  defined  in  (11) 

The  proof  of  this  theorem  is  given  in  Appendix  A.' 

For  nQ  =  0,  condition  (13)  of  Theorem  3  becomes: 

(B°)k(F(A,r)  ]  =  2_kfk(A,r)  >  0  for  every  k  ^  2 

The  above  condition  is  satisfied  if  f2(A,r)  '  0,  since  then 
fk(A,r)  is  a  monotone  increasing  sequence. 

Corollary  2 

In  the  class  of  sequences  satisfying  (12) ,  system  (9) 


has  a  unique  solution  X  =  (xkK 


I 


4 
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0  <  b'k  -  c'  <  .  <bk-c 

—  k  — 

if  f 2 ( A /  r)  =  (1  -  2x ) (r  +2)  -  2X  >  0  (15) 

where 

b  =  g2(r)/f2(X,  r)  =  2 (r  +  3)/((l  -  2X)  (r  +  2)  -  2X) 
and 

X,  r,  b’,  c’,  c  are  as  defined  in  Theorem  3. 


If  the  incoming  traffic  is  Poisson  distributed,  that  is 


then,  condition  (15)  can  be  expressed  in  terms  of  the  traffic 
intensity  X  only,  since  r  =  P(1)/P(0)  =  X. 

Corollary  3 

For  the  Poisson  distribution,  system  (9)  has  a  unique 
solution  X  =  (x^}  satisfying  (9) , 


1-2X  (k~2x)  1  -  xk  -  X+2-2X ( X+3)  (k  2X)  ~  1 


for  X  <  X  =  0.35078 
P 

where  X^  is  the  positive  root  of  the  equation  X+2-2X(X+3) 


=  0, 


For  the  rest  of  the  paper,  we  assume  that  the  input  process 
is  Poisson.  The  following  lemma  gives  an  easily  computable 
recursive  expression  of  the  nc  power  of  the  operator  B,  that 
will  be  very  helpful  in  improving  the  sufficient  condition  given 
in  Corollary  3,  and  in  calculating  the  solution  of  system  (9) 

(if  it  exists) : 


1 
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Lemma  1 


Let  the  sequence  Y  =  {y^ }  be  such  that 

yk  =  (ao  +  b0k)2_k;  k  i  2 


where 


a0,  b0eR 


Then,  for  the  Poisson  distribution 


(Bn).  [Y]  =  l  (a  <i)  +  kb  (i) ) 
k  ^-g  n  n  i 


where 


v0  =  Y'2>  vi+l  =  (1+vi)/2'  a0(0)  =  a0 '  b0(0)  =  b0 


i  (0)  =  -2exp(-X)  V  ((1+Xv.)a  . (i)  +  Xv.b  1  ) 

n  r  i  n-1  l  n-1 


b  (0)  =  -2exp(-X)  y  (a  . (i)  +  b  . (i)) 
n  r  ^_q  n-1  n-1 


a  (i)  =  2exp  (-X  ( 1-v.  .  ))(a  .  (i-1)  +  Xv.  .  b  (i-n)) 

n  l-l  n-1  l-l  n-1 


,  l  <  l  <  n 


b  (i)  =  2exp (-X (1-v.  . ) ) v .  (1  +  v.  .)  b  .(i-1),  1  <  i  <  n 

n  l-l  l-l  l-l  n-1  —  — 


The  proof  of  Lemma  1  is  given  in  Appendix  B. 
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Theorem  4 

For  p(i)  the  Poisson  distribution,  problem  (9)  has  a  unique 
nonnegative  solution  satisfying  (12)  if 

X  <  0.3601 


Proof 

Using  condition  (13)  of  theorem  3,  it  suffices  to  find  some 
nQ  <  <»  such  that: 

n0 

(B  )^[F(A)]  >  0  for  every  k  >_  2  and  A  <  A  =  .3601 

where 

F ( A)  =  { 2_kfk  (X )  } 

fk(A)  =  ( 1-2  A )  k  -  A  (1  +  2  A  )  ;  k  ^  2 

If  A^  <  A2  then  (B  J^lFtA^)]  >_  (B  )^[F(A2)  ]  ,  since,  for 
every  k  _>  2,fk(A)  is  a  monotone  decreasing  function  of  A ,  :  .d 

the  defining  coefficients  (m)  of  the  ope* .  tor  B  ir.  (11)  are 
nonnegative.  Thus,  to  prove  the  theorem,  _t  suffices  to  find 
some  nQ  <  “,  such  that: 

n0 

(B  u )  k  [F  ( A )  ]  >  0  for  every  k  >_  2  (17) 

Let  a^  =  -  A  ( 1  +  2~\)  and  bQ  =  (1  -  2\)  .  Then,  using  Lemma  1 
it  is  not  difficult  to  show  that  (17)  is  true  for  n^  >_  5. 


To  indicate  the  tightness  of  the  above  sufficient  condition, 


n 


we  mention  that  (B  ® ) [F (0 . 3602) ]  remains  negative  for  at  least 


n0  1  75‘ 
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The  next  two  theorems  are  parallel  to  Theorem  4  and  5  in  [12]. 
Theorem  5 

For  p(i)  the  Poisson  distribution,  system  (9)  has  no  nonnegative 
solution  satisfying  (12)  if 

X  >0.363. 

The  proof  of  Theorem  5  can  be  found  in  Appendix  C. 

Theorem  6 

If  system  (9)  has  a  solution  satisfying  (12) ,  then  it 
increases  not  more  rapidly  than  linearly  with  k. 


The  proof  of  Theorem  6  is  omitted. 

We  proceed  with  a  theorem  that  links  the  finite  solution 
of  system  (9)  (if  it  exists)  to  the  solution  {L^ }  of  system  (7) 
for  the  mean  session  length  with  specified  multiplicity  of  the 
algorithm.  This  ^s  necessary ,  since,  formally,  system  (9)  always 
has  the  trivial  solution  x^  =  x^  =  1 ,  x^  =  k  _>  2 .  We  have 

Theorem  7 

If  X  =  (x^ )  is  a  solution  of  system  (9)  satisfying  condition 
(12),  then  =  x^ ,  k  _>  0. 

The  above  theorem  is  parallel  to  Theorem  6  in  [12]  and 
its  proof  is  omitted. 

In  view  of  Theorems  4,  6,  and  7  we  conclude  that  the  average 
session  length  L  given  by  (4)  is  finite  if  X  <  .3601.  But  if  L 

is  finite,  then  the  average  delay  D  is  also  finite,  since  D  <_  L-l. 
Hence,  for  the  Poisson  distribution  the  algorithm  is  stable  if  the 
input  intensity  x  is  less  than  0.3601  packets  per  channel  slot. 


4.  Average  Delax 


In  this  section,  we  calculate  the  mean  session  length  L, 
which  serves  as  an  upper  bound  for  the  average  packet  delay  D 
in  the  stable  region  of  the  algorithm. 

We  first  solve  system  (7)  to  find  the  mean  session  length 
of  specified  multiplicity  for  X  in  the  stable  region.  From 
Theorem  3  we  have 

(An)k[(0)X]  <  Lk  <  (An)k[X(0)]  for  every  n  >_  nQ ;  k  >_  2 

and 

Lk  =  lim(An)R[ (0)X]  =  lim(An)k[x<0)] ?  k  >  2 
n-*-®  n-*-* 

where  A,  nQ,  ^X,  X^  are  as  defined  in  Theorem  3. 

After  simple  calculations,  we  have 

Ak[(0)X]  =  b'k-c '  +  2p(0)b'2~k;  k  >  2 

Ak[X(0)]  =  bk-c  -  2p(0) ( (b(l-2X)-2X)k-X(l+2x)b-2(X+l)2~k 

k  >  2 

where  ^X,  X^,  b,  b' ,  c,  cf  are  as  defined  in  Theorem  3. 

In  view  of  (10)  and  (11)  it  is  not  difficult  to  prove  that 
the  following  are  true: 

<A"M<0,X1  -  b'k-c  ♦  Y  <Bi,k‘,0,Yl  <18> 

K  i=0 


(An).  [X(0))  =  bk-c  -  j  (B1).  [Y(0)] 
K  i  =  0  K 


(19) 
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where 


(0) 


Y 


( (0>ytl: 


,-k 


k  >  2 


and 


a'0  =  2p (0)b* 


(0) 


-  {y 


(0) 


} ; 


(0)  , 
Yk  =  (a 


0 


+  bQk) 2 


-k 


k  >  2 


aQ  =  -2p(0) (X (l+2X)b-2 (A+l) ) 
bQ  =  2p  (0 )  (b  (1-2  A) -2X ) 

The  forms  of  the  lower  and  upper  bound  for  the  solution  L^ 
given  in  (18)  and  (19)  respectively  are  well  suited  for  the 
application  of  Lemma  1.  We  used  Lemma  1  to  calculate  the  powers 
of  the  operator  B  appearing  in  (18)  and  (19).  For  X  <_  0.3  we 
found  that  the  values  of  the  upper  and  lower  bound  coincide  up 
to  the  fourth  decimal  point  within  the  first  fifty  iterations. 

We  should  note  here  that  using  Lemma  1  one  can  calculate  the 
exact  solution  of  system  (7)  for  the  mean  session  length  with 
arbitrary  accuracy  for  any  x  in  the  region  of  stability  of  the 
algorithm. 

In  Table  1,  we  give  the  values  of  the  mean  length  of 
sessions  of  multiplicities  up  to  ten  for  different  values  of  the 
input  intensity. 

The  values  of  L^  given  in  Table  1  were  used  in  (4)  to 
calculate  the  mean  session  length  L.  The  results  are  plotted  in 
Figure  1,  which  gives  the  upper  bound  L-l  for  the  average  packet 
delay  of  the  algorithm 
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5 .  Comparison  of  the  LSBFA  to  Other 
Access  Algorithms 


A  random-access  algorithm  that  can  be  implemented  using 
CNC  binary  feedback  was  first  treated  by  Capetanakis  [4,  5],  and 
by  Tsybakov  and  Michailov  [6].  This  algorithm  uses  the  CTMCRA 
to  resolve  collisions  and  there  are  two  versions  of  it,  a  static 
version  and  a  dynamic  version.  For  the  Poisson  model  and  infinite 
user  population  the  static  algorithm  achieves  a  maximum  stable 
throughput  of  0.346  packets  per  slot  (p.p.s.),  while  the  dynamic 
algorithm  achieves  a  maximum  stable  throughput  of  0.429  packets 
per  slot.  Recently,  Mehravari  and  Berger  [9]  proposed  a 
f irst-come-f irst-served  collision  resolution  algorithm  with  CNC 
binary  feedback,  which  is  stable  for  input  rates  less  than 
0.4422  packets  per  slot.  Also,  Hajek  and  Van  Loon  [15]  have 
recently  shown  that  Aloha-type  retransmission  control  policies, 
that  achieve  a  maximum  stable  throughput  of  e  ^  =  0.3678 
packets  per  slot,  can  be  implemented  on  a  random-access  channel 
using  CNC  binary  feedback. 

All  these  schemes  assume  that  each  user  monitors  the 
feedback  channel  constantly  (for  every  channel  slot  at  all 
times)  even  if  he  has  no  packet  to  send.  This  feedback  requirement 
can  be  reduced  in  the  scheme  proposed  in  [15]  at  the  expense  of 
increased  average  delay. 

In  contrast  to  all  the  schemes  mentioned  above,  the  LSBFA 
eliminates  completely  this  undesirable  and  not  practical 


necessity  by  allowing  continuous  entry  of  new  packets  into  the 
random-access  system.  Furthermore,  the  LSBFA  is  easier  to 
implement  than  all  these  schemes--its  implementation  requires 
only  a  single  counter  possessed  by  each  user. 

In  section  3  we  showed  that  the  LSBFA  is  stable  for  input 
rates  less  than  0.3601  packets  per  slot.  Obviously,  as  the 
level  of  feedback  information  inspected  by  each  user  decreases, 
the  maximum  stable  throughput  decreases  also.  Notice,  however# 
that  the  LSBFA  outerperf orms  the  static  "Tree"  algorithm  [4,  5] 
(0.360  p.p.s.  versus  0.346  p.p.s.),  even  though  the  later  uses 
more  feedback  information. 

The  LSBFA  uses  the  "continuous-entry"  idea  introduced  in 
the  "Stack"  algorithm  [12] .  The  "Stack"  algorithm  uses  ternary 
feedback  (0,1, or  e ) , and  for  the  Poisson-inf inite-population 
model  it  is  stable  for  input  rates  less  than  0.384  packets  per 
slot.  Thus,  in  cases  where  a  central  facility  supplies  the 
feedback  information  by  an  Acknowledgment  scheme,  simplifying 
from  ternary  feedback  to  CNC  binary  feedback  does  not  significantly 
reduce  the  efficiency. 

The  LSBFA  is  easier  to  implement  than  the  "Stack"  algorithm 
because  it  eliminates  the  memory  necessity  in  the  later  algorithm. 
In  the  "Stack"  algorithm  each  blocked  user  at  time  t  has  to 
'Vemember"  the  outcome  of  the  last  nonempty  slot.  The  LSBFA 
requires  no  memory  by  the  users. 

We  now  proceed  to  compare  the  two  algorithms  in  terms  of 
robustness  in  the  presence  of  channel  errors.  The  LSBFA  resolves 
any  collisions  using  the  CTMCRA,  while  the  "Stack"  algorithm 


28 


resolves  any  collisions  using  the  CMTMCRA •  In  [14]  Massey 
showed  that  under  the  more  realistic  situation  where  channel 
noise  can  affect  the  transmissions  on  the  forward  and/or 
feedback  channel,  the  CTMCRA  is  extremely  robust,  while  the 
CMTMCRA  can  suffer  deadlock.  It  is  clear  that  the  "Stack" 
algorithm,  even  though  it  uses  the  CMTMCRA,  does  not  suffer 
deadlock  because  of  the  continuous  entry  of  new  packets  into  the 
random  access  system.  Nevertheless,  we  can  easily  show  that  it 
is  less  robust  than  the  LSBFA.  From  all  the  possible  types  of 
errors  consider  the  one  where  an  "empty"  slot  is  detected  as  a 
"collision"  slot  by  the  blocked  users  because  of  noise  on  the 
channel.  In  the  terminology  of  "Stack"  algorithm,  all  users 
with  a  packet  in  cell  "r"  (r>^l)  of  the  stack  place  their  packets 
in  cell  "r+1".  Thus,  each  time  an  "empty-to-collision"  error 
occurs,  both  cell  "0"  and  cell  "1"  become  empty,  and  the 
algorithm  proceeds  to  resolve  a  nonexisting  collision.  All 
blocked  users  detect  a  sequence  of  empty  slots  but  they  do  not 
move  their  packets  downwards,  since  the  last  nonempty  slot  was 
a  collision  slot  (in  this  case  a  false  collision).  This  deadlock 
situation  lasts  until  some  new  packet (s)  enters  the  system. 

This  happens  with  probability  l-p(0)  independently  at  each  slot. 
Therefore,  the  average  length  of  a  .eadlock  period  is  ld  =  l/(l-p(0)) 
slots.  Under  light  traffic  conditions,  the  average  deadlock 
period  becomes  long  resulting  in  increased  average  packet  delay. 

For  example,  if  p(*)  is  the  Poisson  distribution  with  X  =  0.1,  then 


10.5  slots. 


! 


29 

In  contrast  to  the  "Stack"  algorithm,  the  LSBFA  overcomes 
the  same  error  by  wasting  only  two  slots  independently  of  the 
input  intensity,  as  it  can  be  seen  from  the  rules  of  the  algorithm 
given  in  section  2. 

Using  similar  arguments  like  the  one  used  in  the  "empty-to- 
collision"  type  of  error  case,  we  can  easily  show  that  the  LSBFA 
compared  to  the  "Stack"  algorithm  exhibits  superior  robustness 
in  the  presence  of  channel  errors  of  every  possible  sort. 


k. 


APPENDIX  A 


Proof  of  Theorem  3 

Existence  Consider  the  sequence  X^,  n=0,  1,  2,  .  . 
as  follows: 


,(n)  _  (n)  „ 


=  1 


where 


x^n)  =  (An)klx(0)];  k  >  2 


X{0)  =  {x^0)},  40)  *  bk-c,  k  >  2 


beR  ,  c  =  2bX+l,  and  A  is  the  operator  defined 
For  n=l,  after  simple  calculations,  we  obtain 


XkX)  =  Xk0)  ~  2p (0)  (2  kfkb-2~kgk),  k  >  2 

where 

fkU,  r)  =  (l-2X)k  +  r  ( 1— 2  X )  -  2X 
gk(r)  -  2(k+r+l) 


X  =  l  ip(i) ,  r  =  p(l)/p(0) 
i=0 


From  (A. 1 )  and  (A. 2]  for  k  >  2 ,  we  have 


xkn+1  =  x,[n)  -  2p(0)  (  (Bn)k[F(A,  r)]b  -  (Bn)k[G' 


F  ( x ,  r)  =  {2~kfk(x,  r} 

G(r)  =  {2'kgk(r)  } 

and  the  operator  B  is  as  defined  in  (11) . 


defined 


(A.  1) 


in  (10)  . 


(A. 2) 


r)])  (A. 3) 


where 


31 


If  the  p.m.f.  of  the  input  increment  p(*)  satisfies  condition 
(13)  ,  then 

n0 

(B  )k[F(A,  r)  ]  >  0  for  some  nQ  <  «  and  every  k  >_  2 
Hence,  in  view  of  (A. 3)  we  can  always  find  a  beR+  such  that 

x^n+^  <_  x^n^  for  every  n  >_  n^  and  k. 

For  a  given  p(-),  let 

n0  nn 

b  =  max(  (B  u)k[G{r)  ]/(B  u)k[F(X,  r)  ])  ,  k>2  (A. 4) 

Jc 

Clearly,  if  b  is  chosen  as  in  (A. 4),  then  for  every  k  >  2, 


(n0 ) 

the  values  ^ 


(nQ+1 )  (nQ+2) 


form  a  nonincreasing 


sequence- 

If  A  <  1/2,  then  *k  ;  is  nonnegative  for  every  k  and  n, 

since  with  b  chosen  as  in  (A. 4)  *k°^  is  nonnegative  for  every  k 

and  the  operator  A  is  nonnegative  (qk  ( • )  _>  0)  .  Note  that  if 

A  _>  1/2,  then  there  is  no  nn  for  which  (13)  is  true. 

(nQ)  (nQ  +1 )  U(nQ+2) 

Since  xk  '  '  Xk  ,  is  a  nonnegative 

nonincreasing  sequence  the  following  limit  exists  for  every  k  >  2 


x  =  lim  x 
k  k 

n-*oo 


(A. 5) 


Clearly,  if  (13)  is  true,  then  0  <_  x^n  <  00  for  every  n  n^  and 
2  <_  k  <  <*>.  Moreover  for  every  n  >  n^  we  have 


,(n+i)  a  r  _  x(n) 

V  -  i  +  2  l  qk  (m)  xm  <  x, 


Thus,  for  every  n  and  2  <  k  <  ®  the  nonnegative  series 


2  <Jkln,)1‘mr 
m=u 
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is  absolutely  convergent  (satisfies  condition  (12)).  Therefore, 
for  arbitrary  e  >  0  and  n  there  exists  a  Kq  ,  such  that  foi 
arbitrary  K  >kQ  we  have: 


m=K  qk{ir) 


x  <n> 
m 


<  e/10, 


k  <  K 


(A. 6) 


In  view  of  (A. 5)  and  for  the  given  e  and  K  we  have 

|x^n)  —  Xj,,  J  =  xjjn)  -  xk  <  e/5,  for  all  n  >  N  _>  nQ  and  k  <  K 


Then,  for  arbitrary  kQ  <  »  we  obtain 

K-l 

l*k  "  Ak  [X]  I  i  lxv  “  2  l  (m)xm  "  Xl  +  |  1 
K0  K0  K0  m=0  K0  m  5 


(A. 7) 


(n+1) 


K-l 


-  2  l  qk  xmn)  “  +  b  -  lx 

m=0  K0  m  D 


^n+1)  -  A.  (X(n)]|  +  e 
0  K0 


Hence  X  =  (xk)  is  a  solution  of  system  (9) ,  since  the  above 
inequality  is  true  for  arbitrarily  small  e  >  0. 

We  proceed  now  with  the  lower  bound  sequence  ^n^X,  n=0,  1,  2,  . 
which  is  defined  as  follows: 


where 


<n>x  =  <n>x  =  i 
x0  X1 

(n,xk  =  (An)k[(0)X],  k  >  2 


(0>  X  =  {  (0)x  },  (0)x  =  b'k-c '  ,  k  >  2 

k  k 

b ’  =  2/  (1-2X)  ,  c'  =  2b' X  +  1 


for  n=l  after  simple  calculations,  we  obtain 
(1)x.  =  (0)  x.  -  2p(0)2_kb' 


(A. 8) 
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In  view  of  the  fact  that  the  operator  A  is  nonnegative 
(q^(*)  ±  0) ,  it  follows  from  (A. 8)  that  for  \  <  1/2  and  f^r 
every  k  _>  2,  the  values  ^x^,  ^x^,  ^x^,  .  .  .  form  a 
nondecreasing  sequence.  This  sequence  is  also  bounded  from 
above  by  the  sequence  X^.  Indeed,  if  condition  (13)  is  true, 
then  b'  <  b,  where  b  is  as  defined  in  (A.  4)  .  (If  b'  >_  b  then  it 
follows  from  (A. 2)  that  Xj^_>  for  every  k  _>  0  and  consequently 

condition  (13)  cannot  be  true.)  Thus,  ^x^  =  b'k-c'  <  bk-c  =  x^^ 
for  every  k  _>  2  and  since  the  operator  A  is  nonnegative  it 
follows  that 


Since  ^ 


(Dx 
'  \ 


for  every  n  and  k  >_  2  . 

(2) 

x^,  .  .  .  is  a  nondecreasing  sequence, 

(n) 


which  is  bounded  from  above  by  the  sequence  x  ,  the  following 
limit  exists  for  every  k  _>  2: 


1  (n) 

\  •  i1"  \  i  \ 


(A. 9) 


Following  the  same  steps  as  in  the  case  of  X  =  { x^ }  we  can 

r 

prove  that  the  sequence  X'  =  { x^ }  is  a  solution  of  system  (9) 
satisfying  (12) . 

Uniqueness  To  prove  the  uniqueness  of  a  solution  we  assume 
that  system  (9)  has  two  distinct  solutions  satisfying  (12)  and 
we  end  up  with  a  contradiction.  The  proof  parallels  the  one  given 
in  [12,  p.  236]  and  it  is  omitted. 


From  the  uniqueness  of  t*~°  solution,  we 

- lim(n)  *k  - lim  4n)  -  \ 

n-*-®  n->® 


have 
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APPENDIX  B 


Proof  of  Lemma  1 

We  prove  lemma  1  by  induction.  For  n=l , 


we  have 


I  VY1 


-  I  qk(m)  (a0  +  b0m)vj 


=  a. 


l 

m=0 


m 


qk  (m)  vQ 


+  b. 


L 

m=0 


m 


qk(m)mv0 


'  <a0qk(0) 


+  2(a0+b0)qk(1 


(B.l) 


where  qk (  • )  is  as  defined  in  (8)  and  Vq  =  0.5.  In  particular, 
qk(0)  =  exp  (~  A)  v^ ,  and  qk<l)  =  exp  (- x)  (k+  x)  v* .  After  simple 
calculations,  we  obtain 


00  1  4- V 

l  q  (m)vm=  exp  (- >.  (1-v) )  (— j-) 
m=0  K 


(B .  2 ) 


£  qk(m)mvm=  exp(-X  (1-v)  )  k+\v)(— j-*)  (B 

m=0 

Substitution  of  qk(0),  q^d),  (B.2),  and  (B.3)  with  v=vQ  into 

(B.l)  gives 

Br[Y]  =  ax  (0)  Vq  +  ax  (1)  v^  +  b^OkvJ  +  bx(l)kv^  (B 


where 

vx  =  (l+v0)/2, 

3l(0)  =  -2exp(-x)(aQ  +  *v0(a0+b0))' 
b1(0)  =  -2exp(-x)v0(aQ+b0) , 

a^l)  =  2exp(-x  (l-vQ) )  (a0  +  xvQb0)  , 
b:(l)  =  2exp (-X ( 1-vq^  vq (l+vQ) _1bQ 
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i 


i 


Equation  (B.4)  proves  the  lemma  for  n=l. 

Next  we  assume  that  (16)  is  true  for  n=j  and  we  prove  that 
it  is  true  for  n=j+l  as  well.  We  have 

(B3).  [Y]  =  l  (a.  (i)  +  kb  (i))v* 

K  i=0  3  3  1 


(B3  +  1)  [Y]  =  B,  [B3  [Y]  ]  =2  l  q.  (ip.)  \  (a  (i)  +  mb.  (i))v® 

k  k  m=2  k  i=0  3  3 


=  2(  [  a.(i)  l  q.  (m)  v 

i=0  3  m=0  J'  i=0  J  m=0 


+  j 

i  .  L, 


m 


k . i  ■  ,i„bj(i'_L  V"”"vi 


-  (qv(0) 


i=0 


a  j  (i)  +  qk  (i) 


l  (a  .  (i)  +  b  .  (i)  v  ) 
i=0  3  3 


(B.  5) 


Substitution  of  qk(0),  qk(l),  (B.2)  and  (B.3)  with  v=vi  into 
(B. 5)  gives 


(Bj+1).  [Y]  =  'l  (a  (m)  +  kb  (m)) 


j  +  1 


1=0 


j  +  1 


j+1  '  ' '  m 


(B .  6 ) 


where  a ...  (it)  ,  b .  , .  (ir)  ,  m=0 ,  1,  .  .  .,  j  +  1,  are  as  given  in  (16). 
3  +  x  3  i 

Equation  (B.6)  proves  the  lemma  for  n=j+l. 


4 


i 
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APPENDIX  C 

Proof  of  Theorem  5 

Consider  the  following  truncated  system  of  equations 

*k  =  1  +  2J„  V">V 

m=0 

or  in  matrix  notation 


where 


-  Q  )x  = 
n  -n 


c 

— n 


(C.l) 


-n  “  (x2'  X3'  ‘  r  Xn+1)T 

-n  =  (c2 '  C3 '  •  *  "  cn+l}  '  ck  =1+  qk(0)  +  qk(1)' 
k  =  2 ,  3 ,  .  .  . ,  n+1 

Q  =  (q. .)  is  a  (nxn)  nonneaative,  irreducible,  square 
n  ^i] 

matrix  with  q^  =  q^(j),  2  <_  i  £  n+1,  2  £  j  <_  n+1, 
and  I  is  the  (nxn)  unit  matrix 


The  following  theorem  gives  the  conditions  that  ensure 
positivity  (x  >  0)  of  solutions  to  the  equation  system  (C.l). 


Theorem  [16,  Th .  2.1] 

A  necessary  and  sufficient  condition  for  a  solution 
xn(xn  >  0,  /  0)  to  equations  (C.l)  to  exist  for  any  £n  0,  /  0 
is  that  r^  <  1/2;  where  rn  is  the  Perron-Frobenius  eigenvalue 
of  the  nonnegative,  irreducible  matrix  Qn .  In  this  case,  there  is 
only  one  solution  xn,  which  is  strictly  positive  and  given  by 


x 

— n 


c 

-n 


■SCTra 


3? 


It  is  known  from  the  theory  of  the  nonnegative  matrices  (see 
[16])  that 


'n+1 


n 


and  lim  r  =  r 
n 

n-*-« 


(C .  2) 


where  r  is  the  Perron-Frobenius  eigenvalue  (the  reciprocal 
"convergence  norm")  of  the  infinite  dimensional  matrix  Q;  where 
Q=lim  Qn  as  n-*-°°  (Qn  is  the  (nxn)  northwest  corner  truncation  of 
Q)  • 


From  the  above  theorem  and  (C.2)  it  follows  that  if  for 

some  nn  we  have  r  >  1/2  then  for  all  n  >  n„  system  (C.l)  has 
u  nQ  -  u 

no  nonnegative  solution.  In  this  case,  equation  system  (9)  has 
no  nonegative  solution  satisfying  (12)  ,  since  it  is  obtained 
from  equation  system  (C.l)  as  n->-°°.  To  calculate  the  Perron- 
Frobenius  eigenvalue  rn  we  made  use  of  the  following  lemma: 


Lemma  [ 17 ] 

If  s(T)  and  S (T)  are  the  minimal  and  maximal  row  sums  of  a 
square  nonnegative,  irreducible  matrix  T  with  Perron-Frobenius 
eigenvalue  r,  then 


s (T)  <  (s(T2))1/2  <  (s(T4))1/4  <  .  .  .  <  r  <  .  .  .  < 


<(S(T4))1/4  <  (S(T2))1//2  <  S(T) 


Calculations  made  for  n^  =  10  show  that  r^g  >  1/2  for  X  =  0.363. 
Thus,  for  \  >_  0.363  equation  system  (9)  has  no  nonnegative 
solution  satisfying  (12)  . 
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Table  1.  Lk(X)  :  Mean  length  of  a  session  of  multiplicity  k 
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members,  a  majority  of  whom  conduct  research  in  addition  to  teaching. 

Research  is  an  integral  part  of  the  educational  program  and  interests  parallel  academic  specialties.  These 
range  from  the  classical  engineering  departments  of  Chemical,  Civil,  Electrical,  and  Mechanical  and 
Aerospace  to  departments  of  Biomedical  Engineering,  Engineering  Science  and  Systems,  Materials 
Science,  Nuclear  Engineering  and  Engineering  Physics,  and  Applied  Mathematics  and  Computer  Science. 
In  addition  to  these  departments,  there  are  interdepartmental  groups  in  the  areas  of  Automatic  Controls  and 
Applied  Mechanics.  All  departments  offer  the  doctorate;  the  Biomedical  and  Materials  Science 
Departments  grant  only  graduate  degrees. 

The  School  of  Engineering  and  Applied  Science  is  an  integral  part  of  the  University  (approximately  1 ,530 
full-time  faculty  with  a  total  enrollment  of  about  1 6,000  full-time  students),  which  also  has  professional 
schools  of  Architecture,  Law,  Medicine,  Commerce,  Business  Administration,  and  Education.  In  addition, 
the  College  of  Arts  and  Sciences  houses  departments  of  Mathematics,  Physics,  Chemistry  and  others 
relevant  to  the  engineering  research  program.  This  University  community  provides  opportunities  for 
interdisciplinary  work  in  pursuit  of  the  basic  goals  of  education,  research,  and  public  service. 


